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INTRODUCTION 

Assume  an  explosive  source  located  at  it,  with  coordinates  X and  Y, 
occurring  at  time  T.  The  resulting  wave  may  be  described  by  rays  which  are 
normal  to  the  wave  front  and  which  travel  with  velocity  cn  + w,  where  c is  the 
mean  sound  speed,  n is  the  unit  vector  in  the  direction  of  propagation  relative 
to  the  air,  and  w is  the  wind  velocity,  with  x-  and  y-components  u and  v.  The 
wave  arrives  at  a microphone  located  at  r^  = (x^,  y ) at  time  t^.  The 

following  vector  equality  holds: 

R + (cini  + wi)  (ti  - T)  = ?1,  (1) 

where  c^  and  w^  are  mean  speeds  over  the  path  from  the  source  to  M^.  When  Eq. 
(1)  is  solved  for  cn^(t^-T),  and  resolved  into  component  form, 

Ci(ti  ' T)2  = -X-xi  + ui(ti  " T)l2 

+ [Y  - y±  + vi(ti  - T)]2.  (2) 

Formally,  one  has 

fi(X,Y,T)  = 0,  i = 1,2,. ..n,  (3) 

if  there  are  n microphones. 

The  system  of  non-linear  equations  represented  by  Equation  (3)  requires 
three  microphones  for  a solution  of  X,  Y and  T.  More  than  three  microphones  form 
an  over-determined  system.  To  completely  locate  the  source  requires  accurate 
knowledge  of  c^,  w^,  r^,  and  t^. 

Devoting  attention  to  the  problem  of  determining  t^  alone,  without  the 
consideration  of  the  complete  problem  could  be  merely  an  academic  exercise. 





First,  it  is  obvious  that  the  error  in  determining  the  source  position  is 
dependent  on  all  the  quantities  c^,  w^,  r ^ , and  t^.  A surveying  error  of  one 
foot  toward  the  source  is  equivalent  to  one  msec  time  error.  The  quantities 
c^  and  w.  must  be  considered  means  over  the  individual  ray  paths,  which  do  not 
necessarily  follow  the  surface.  Second,  the  problem  of  what  is  the  "best" 
solution  to  an  over-determined  system  of  equations  is  not  clear.  The  present 
method  appears  to  weight  the  end  microphones  by  one-half.  Should  a weighted 
least-squares  method  be  used,  or  perhaps  a "robust"  method  which  is  insensitive 
to  large  errors  in  individual  microphones?  An  empirical  test  of  the  "best" 
determination  of  t . is  not  independent  of  the  algorithm  used  in  solving  Equation 
(3).  Third,  the  nature  of  the  signal  is  changed  by  the  propagation  through  the 
air  and  by  the  detection  and  recording  system.  This  fact  severely  affects  the 
determination  of  t.. 

i 

Thus  certain  areas  of  the  entire  sound  ranging  problem  — those  areas 
which  appear  to  directly  affect  the  determination  of  arrival  times  — will  be 
investigated.  The  total  problem  is  shown  as  a block  diagram  in  Figure  1.  The 
source  is  characterized  by  its  position  X and  Y,  time  T and  energy  E.  This 
produces  some  disturbance  in  time  f(t).  This  disturbance  is  propagated  through 
the  atmosphere,  which  acts  as  a non-linear  filter,  producing  a different  signal 
h(t).  This  filter  may  be  characterized  by  its  sound  speed  c,  a function  of 
time,  position,  and  height  and  by  the  ratio  of  the  specific  acoustic  impedance 
of  the  surface  to  that  of  air,  Z.  The  microphone  detects  this  signal  and  what- 
ever noise  is  available.  The  microphone,  the  amplifier,  and  the  recorder  also 
act  as  a filter,  with  response  g(w),  changing  the  signal  to  q(t).  It  is  this 
signal  which  provides  the  information  from  which  the  arrival  time  at  each  micro- 
phone must  be  determined.  The  arrival  times,  together  with  surveying  and 


Figure  1.  Block  diagram  of  sound  ranging  problem.  The  dashed 
lines  indicate  the  theoretical  connections  by  means  of  a model 
At  present,  there  are  no  connections  from  source  or  detection. 
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meteorological  data,  are  analyzed  according  to  a model  which  is  supposed  to 
represent  the  propagation  delay  time.  The  result  is  the  position  and  perhaps 
the  strength  of  the  source. 

Section  I will  deal  with  the  nature  of  the  source  signal  during  the  early 
part  of  the  propagation,  when  it  must  be  treated  as  a shock.  Section  II  con- 
siders the  propagation  problems  arising  from  vertical  gradients  in  temperature 
and  wind,  and  the  ever-present  boundary  between  the  earth  and  air.  The  effect 
of  the  microphone  response  on  the  signal  will  be  treated  in  Section  III. 

Finally,  Section  IV  will  contain  recommendations  for  future  investigations. 

Some  attempts  at  digital  analysis  and  semi-automation  in  determining  arrival 
times  will  be  described  in  another  paper,  although  some  general  comments  on  this 
subject  are  included  in  the  recommendations. 


I 

r 
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I.  THE  PROPAGATION  OF  AN  EXPLOSIVE  SHOCK  WAVE 
The  propagation  of  explosive  shocks  in  the  atmosphere  has  received 
extensive  investigation  during  the  last  30  years.3  Both  theoretical  calcula- 
tions and  experimental  measurements  over  wide  ranges  of  explosive  energies 
agree  that,  after  a certain  distance,  the  disturbance  is  described  by  the  pres- 
sure vs  length  signal  shown  in  Figure  2.  Both  the  magnitude  of  the  pressure 
rise  P and  the  length  of  the  positive  pressure  pulse  L are  functions  of  the 
dimensionless  variable  £ = r/d,  where  r is  the  distance  from  the  source  and  d 
is  a characteristic  length  given  by 

d = (E/P  )1/3,  (4) 

0 

where  E is  the  energy  of  the  source  and  p^  is  the  ambient  atmospheric  pressure. 

Less  is  known  about  the  negative  pressure  region.  The  magnitude  and 

2 

length  described  in  Figure  2 are  taken  from  Brode  . Sometimes  secondary  reflect- 
ed pulses  ride  along  this  region.  However,  the  initial  pressure  surge  and  the 
length  of  this  pulse  are  certainly  the  important  parameters  of  the  signal  in 
determining  arrival  time.  Surely  if  this  signal  were  not  distorted  by  the 
atmosphere,  ground,  and  microphone,  the.  determination  of  arrival  times  would  be 
simplified . 

The  characteristic  length  given  by  Equation  (4)  is  proportional  to  the 

cube  root  of  the  energy  of  the  charge  propelling  the  projectile  less  the  energy 

of  the  projectile  and  recoil,  thus  approximately  proportional  to  the  diameter 

of  the  bore.  However,  for  artillery,  various  lengths  of  charge  are  available. 

3 

Using  values  reported  by  Wurtele  and  Roe  , the  energy  for  the  No.  7 charge  for 

g 

an  8-in  bore  is  1.22  x 10  J.  For  sea-level,  this  results  in  d = 10.6  m.  Other 
charges  and  bores  of  US  artillery  have  energies  which  may  be  approximately  by 


i 
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scaling  to  the  kinetic  energy  of  the  projectile,  given  by  the  weight  and  muzzle 
velocity.  Table  1 contains  the  resulting  characteristic  lengths  for  charges  N 
No.  1 through  No.  7 for  8-in,  155-mm  and  105-mm  artillery.  It  may  be  noted  that 
there  is  overlap,  depending  on  the  charge  used,  between  light  and  medium  and 
between  medium  and  heavy  artillery. 


Table  1.  Characteristic  lengths  for  various  charges  and  diameters.  (Length  in 
meters . ) 


Diam. 

Charge  #1 

Charge  #3 

Charge  #5 

Charge  #7 

105  mm 

2.9 

3.3 

3.9 

5.2 

155  mm 

4.1 

4.9 

6.0 

8.0 

8 in 

5.9 

6.8 

8.4 

10.6 

At  the  range  of  1 km.  Table  1 yields  a dimensionless  parameter  4 which 
varies  between  94  for  the  most  powerful  8-in  charge  and  350  for  the  least  power- 
ful 105-mm  charge.  The  largest  value  expected  would  be  for  No.  7,  105-mm  at 
maximum  range  (11.5  km),  which  is  r,  = 2200.  Thus  the  range  of  4 of  interest  in 
sound  ranging  is  from  100  to  2000.  For  these  values  of  4,  the  shock  is  weak  and 
the  compilation  by  Baker  is  incomplete;  his  values  for  L beyond  4 = 60  being  an 
extrapolation.  Fortunately  the  shock  is  weak,  allowing  the  calculation  of  wave- 
shapes and  arrival  times. 

The  detail  of  this  calculation  are  relegated  to  Appendix  A.  Due  to  the 
quasi-adiabatic  nature  of  the  propagation,  the  high-pressui e shock  front  travels 
faster  than  the  ambient  pressure  point  on  the  wave.  The  rate  at  which  L increases 
with  distance  of  propagation  is  given  approximately  by 


8 


I 


where  Y is  the  ratio  of  specific  heats  (1.40  for  air)  and  P is  the  peak  over- 
pressure. The  increase  in  entropy  across  the  weak  shock  causes  the  total  energy 
of  the  wave  to  decrease  with  the  distance  of  propagation  according  to 

l dE  = y+i  P_  1_  ,,, 

E dr  “ 4Y  p L ’ ' 

o 

where  E is  the  energy  of  the  wave,  given  by 

2 2 

r = 4TTr~P~L 


here  p is  the  air  density  and  c is  the  sound  speed. 

The  substitution  of  Equation  (5)  into  Equation  (6)  and  the  integration 


of  the  result  yields 


EL  = constant. 


Equation  (7)  and  (8)  may  be  solved  for  P as  a function  of  r and  L, 


P = P (r  /r)  (L  /L), 
o o o 


where  the  zero  subscripts  refer  to  some  reference  distance  from  the  source,  r 

c 

Finally,  Equation  (5)  combined  with  Equation  (9)  may  be  integrated  to  obtain: 

L = /T[L  +^—  r ln(r/r  )]^2.  (10) 

o o 4y  Pq  o o 

Using  Baker's  tabulated  values  for  z,  = 60,  (P  / p = 2.48  x 10  \ L = 0.8561, 
° * o o o 

r = 60d),  one  has 
o 


L(c)  = d [0. 733  + 0.109  In  (c/60)] 


Several  values  of  this  function  and  others  of  interest  in  sound  ranging 
are  tabulated  in  Table  2.  The  relative  overpressures  P/po  are  determined  from 
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Equation  (9).  The  next  column  contains  the  instantaneous  speed  of  the  wave 
(V/c)  as  a function  of  4.  This  is  calculated  from 

V = c(l  + ^ — ) , (12) 

4Y  Pq 

derived  in  Appendix  A.  One  notes  that  this  speed  is  very  close  to  the  speed  of 
sound . 


The  important  speed  in  the  sound  ranging  equations  is,  however,  not  the 
instantaneous  speed,  but  the  mean  speed  over  the  path.  Due  to  the  strong 
shock  propagation  close  to  the  source,  the  mean  speed  is  larger  than  the  instan- 
taneous speed.  The  last  column  of  Table  2 lists  this  mean  speed  as  a function 
of  4.  For  4 = 60,  Baker  gives  for  the  arrival  time  t = 58.8c,  or  <V>/c  = 1.020. 
The  mean  speeds  in  Table  2 are  based  on  this  plus  the  change  in  arrival  time  due 
to  AL  between  4 values.  The  tabulated  values  are  given  quite  accurately  by  the 
empirical  fit 

<V>/c  = 1 + 1.4/e  - 15/e2  (13) 

over  the  range  e = 100  to  C = 2000.  It  may  be  noted  that  these  values  are  not 
negligible  for  heavy  artillery  at  short  ranges,  varying  between  a 1.25%  correc- 
tion to  a 0.39%  correction  for  heavy  to  light  artillery  at  1 km.  This  is  equi- 
valent to  temperature  corrections  of  1.9°C  and  0.2°C  respectively. 


Table  2.  Dimensionless  variables  for  the  propagation  of  explosive  shocks.  The 

1/3 

characteristic  length  is  d = (E/p^)  . The  characteristic  speed  is 

the  sound  speed  c. 


4 = r/l 

L/d 

P/P 

o 

v/c 

<v>/c 

100 

0.893 

1.43xlO~3 

1.00061 

1.0125 

200 

0.941 

6.77xl0_4 

1.00029 

1.0065 

500 

1.002 

2.54 

1.00011 

1.0027 

1000 

1.045 

1.22 

1.00005 

1.0014 

2000 

1.086 

5.86xl0_5 

1.00003 

1.0007 

10 


In  principle,  a faithfully  received  signal,  together  with  a known  range, 
would  allow  the  determination  of  d and  hence  the  energy  of  the  source.  In 
practice,  the  vagaries  of  propagation  close  to  the  earth's  surface  make  this 
determination  exceedingly  difficult.  The  most  useful  information  is  probably 
contained  in  the  period  or  frequency  content,  rather  than  the  pulse  height. 
Assuming  a mean  value  of  L/d  to  be  1,  the  period  of  the  positive  pulse  will  be 
approximately  d/c.  The  characteristic  lengths  in  Table  1 show  that  heavy  artil- 
lery will  have  periods  ranging  from  17  to  31  msec,  medium  artillery  will  range 
from  12  to  24  msec,  and  light  artillery  will  range  from  9 to  15  msec  in  period 

• **  ■ * • • m • . 

in  the  order  of  12  msec. 

The  frequency  spectrum  of  a signal  such  as  shown  in  Figure  2 is  obtained 
from  the  Fourier  transform  of  the  pulse. 


P (to) 


iu)t 

e 


(14) 


where  to  is  the  angular  frequency  2irf.  For  a pulse  approximating  Figure  2, 


P(t) 


t<0 


P (1-t/r), 

O 

o 


0<  t < T , 
t>T 


(15) 


P (to) 

= __1_ 

p I 

CUT 

o 

(16) 


This  relative  spectrum  is  shown  in  Figure  3.  The  high  frequency  spectrum  falls 
off  at  6 db  per  octave,  being  down  3 db  at  wx  = 5.  This  3-db  point  corresponds 
to  a frequency  of  0.8/x,  where  x is  the  over-pressure  period.  This  frequency 
ranges  from  26  to  46  Hz  for  heavy,  34  to  66  Hz  for  medium,  and  52  to  94  Hz  for 


light  artillery. 


Figure  3.  The  magnitude  of  the  pressure  spectrum  of  a pulse  such  as 
shown  in  Figure  2.  The  period  of  the  pulse  is  T - L/c . The  high 
frequency  rolls  off  at  6 db/octave  with  the  -3  db  point  at  about 
ofT  = 3. 
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The  above  theoretical  predictions  for  the  wave  shapes  due  to  various 
artillery  assume  spherical  propagation  through  a homogeneous,  non-absorbing 
medium  with  no  boundaries.  Perhaps  the  two  most  dangerous  assumptions  involve 
the  neglect  of  absorption  and  the  neglect  of  boundary  influence.  Both  of  these 
neglected  effects  tend  to  filter  out  the  high  frequencies  of  the  wave,  although 
the  fundamental  over-pressure  period  ought  not  be  effected  much.  These  effects 
then  primarily  increase  the  rise  time  of  the  initial  pressure  pulse,  although 
the  finite  amplitude  effect  considered  here  will  tend  to  reshape  the  shock. 

Since  the  neglected  effects  do  not  change  the  period  of  the  pulse,  the 
fundamental  frequency  of  the  period  would  perhaps  be  more  indicative  than  the 
3-db  points  of  the  unabsorbed  pulse.  These  frequencies  range  from  16  to  20  Hz 
for  heavy,  21  to  41  Hz  for  medium,  and  33  to  59  Hz  for  light  artillery.  Section 
III  will  show  that  the  response  of  the  T-23  microphone  is  3db  down  at  about  25  Hz. 
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II.  ATMOSPHERIC  PROPAGATION  OF  THE  PULSE 
The  sound  ranging  equation,  Eq.  (2)  is  a two-dimensional  equation, 
assuming  circular  propagation  in  a plane.  This  plane  is,  of  course,  the  surface 
of  the  earth  where  both  the  source  and  receivers  are  located.  Disregarding  that 
this  surface  is  not  often  a plane  to  within  the  wavelengths  involved,  and  assum- 
ing a homogeneous  surface  temperature  and  wind,  there  still  exists  the  problem 
of  temperature  and  wind  variations  with  altitude.  It  is  well  known  that  these 
variations  cause  waves  to  be  refracted  — sometimes  upward  resulting  in  the 
microphone  being  in  a shadow  zone,  sometimes  downward  resulting  in  direct  propa- 
gation along  a curved  path.  (See  Figure  4.)  Thus  two  cases  exist:  If  the 
phase  velocity  of  sound  decreases  with  altitude,  only  diffracted  energy  reaches 
the  microphone,  traveling  with  the  speed  determined  by  the  surface  temperature 
and  wind.  As  is  true  for  all  diffracted  waves,  the  low  frequency  components  are 
diffracted  with  greater  efficiency  than  the  high  frequency  components.  The 
second  case  of  refracted  propagation  occurs  if  the  phase  velocity  increases  with 
altitude.  Due  to  Fermat's  principle,  the  curved  path  results  in  a travel  time 
which  is  less  than  the  path  along  the  surface.  This  results  in  a higher  mean 
velocity  than  that  predicted  by  use  of  the  surface  temperature  and  wind.  In 
addition,  reflections  from  the  ground  make  possible  multiple  paths,  with  the 
corresponding  superposition  of  signals. 

The  diffracted  case  is  mathematically  difficult  and  is  also  beset  with 
experimental  difficulties  involving  the  determination  of  the  specific  acoustic 

impedance  of  the  surface,  not  to  mention  the  large  number  of  different  geologi- 

3 

cal  and  biological  surfaces  imaginable.  Wurtele  and  Roe  , using  Doak's  solu- 
4 

tion  , have  determined  the  effect  of  the  boundary  on  the  shape  of  pulses  described 
by  Heaviside  step  functions  and  "N"  shaped  waves.  Their  results  indicate  the 


Figure  4.  Propagation  for  (A)  gradient  in  phase  velocity  is  negative,  and 
(B)  gradient  in  phase  velocity  is  positive.  Case  A creates  a shadow  zone 
along  the  surface  however  sound  is  diffracted  (with  loss  of  high  frequency 
components)  along  the  surface.  Case  B results  in  the  refracted  ra'  arriv- 
ing before  a surface  ray  would  arrive. 
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rather  drastic  filtering  effects  on  pulses  diffracted  over  boundaries  which 
have  specific  acoustic  impedances  that  are  not  much  larger  than  that  of  air. 
Figure  5 (adapted  from  their  Figure  A. 5)  shows  the  resulting  shape  of  an  "N" 
shaped  wave  of  29  msec  duration  after  being  diffracted  one  km  along  surfaces  of 
various  specific  acoustic  impedances.  In  the  figure,  Z is  the  ratio  of  the 
surface  impedance  to  that  of  air,  pc. 

Although  the  signal  is  fairly  faithfully  reproduced  for  Z = 100,  lesser 
values  of  Z show  increasingly  severe  loss  of  high  frequency  components.  The 
boundary  thus  acts  as  a low-pass  filter,  reducing  the  amplitude  of  the  signal 
and  increasing  the  uncertainty  in  determining  the  arrival  time,  although  the 
original  pulse  length  is  preserved.  The  rise  time  increases  with  distance,  as 
shown  in  Fig.  6,  where  a unit  step  function  is  shown  after  being  propagated 
according  to  Wurtele  and  Roe  for  distances  of  1,  2,  5,  and  10  km  along  a surface 
with  Z = 10.  Consideration  of  the  fact  that  artillery  pulses  are  not  longer 
than  30  msec  leads  to  the  conclusion  that  arrival  times  are  increasingly  diffi- 
cult to  measure  for  longer  distances. 

Even  though  diffracted  signals  may  (dependent  upon  Z and  range)  make 
source  strength  determination  and  arrival  time  difficult  to  measure,  the  propa- 
gation along  the  surface  simplifies  the  determination  of  mean  sound  speed.  The 
proper  values  to  use  when  the  phase  velocity  decreases  with  height  are  simply 
the  surface  values  of  temperature  and  wind  velocity.  There  appears  to  be  no 
need  to  include  upper-level  values  in  some  weighting  scheme. 

In  cases  where  the  phase  velocity  increases  with  height,  the  situation 
is  reversed.  The  refracted  waves  ought  to  faithfully  represent  the  pulse  in- 
vestigated in  Section  I,  but  the  mean  sound  speed  and  wind  velocity  must  include 
the  properly  weighted  upper-level  meteorological  data.  These  mean  values  could 


Figure  b.  Rise  times  for  a Heaviside  step  function  diffracted  along 
boundary  of  Z = 10  for  various  distances.  (Adapted  from  Wurtele  and 
Roe,  1977). 
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be  determined  by  a ray-tracing  program,  but  that  would  probably  require  too 
much  computer  power.  An  approximation  to  the  correct  weighting  may  be  obtained 
for  simple  atmospheric  models.  The  simplest  model  is  one  in  which  it  is  assumed 
that  both  the  sound  speed  and  wind  vectors  have  linear  gradients  with  respect 
to  height. 

The  analytic  solution  for  such  a model  is  contained  in  Appendix  B.  If 
a Taylor  expansion  for  the  sound  speed  c as  a function  of  height  z is  truncated 
at  the  first-degree  term, 


c (z ) = co  (1  + az ) , 


(17) 


where 


1 dc 

a = — -j-  , 

c dz 
o 

and  c is  the  sound  speed  at  the  surface.  Likewise, 
resolved  into  a component  tangent  to  the  path  of  the 
mal  to  the  path  (v),  may  be  expanded  to  yield: 


(18) 

the  wind  vector,  when 

ray  (u)  and  a component  nor- 


where 


u (z)  = u + c bz , 
o o 

(19) 

v(z)  = v + c dz, 
o o 


1 du  , _ 1 dv 

c dz  ’ c dz 

o o 


Under  the  assumption  that  higher-order  terms  are  negligible,  it  is  shown  in 
Appendix  B that  the  mean  speeds  for  a refracted  ray  are  given  by  the  values  c, 
u,  and  v evaluated  at  one-third  of  the  maximum  height  of  the  ray,  and  that  the 
maximum  height  obtained  by  the  ray  is  given  by 


Z - J (a  + b)  R2, 


(20) 
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where  a and  b are  the  gradients  defined  above,  and  R is  the  range  of  the  ray 
(the  distance  between  source  and  receiver).  It  is  clear  from  Equation  (20)  that 
Z is  positive  only  if  a + b > 0.  Otherwise  the  ray  is  diffracted  along  the 
boundary  rather  than  refracted  along  a curved  path. 


Table  3.  Meteorological  variables 

for  first  example. 

Parameter 

Surface 

800  m-level 

Virtual  temp. 

280  K 

275  K 

East  wind 

0 

5 m/sec 

North  wind 

0 

2 m/sec 

As  an  example,  suppose  the  virtual  temperature  and  wind  components  are 

1/9 

given  by  Table  3.  One  has  C = 20.05  (280)  “ = 335.5  m/sec,  conn  = 332.5  m/sec, 

O oUU 

a = -1.12  x 10  ^ m \ b = -1.86  x 10  ^ m (along  the  x-axis),  and  d = +0.75  x 
10  ^ m 1.  Propagation  in  most  directions  is  along  the  surface,  with  c = 

335.5  m/sec.  However,  propagation  in  the  negative  x-axis  (toward  the  west, 
down  wind  from  the  east  wind)  is  refracted.  Here  a + b = 0.74  x 10  m ^ , yield- 

ing a maximum  height  of  only  one  meter  for  a range  of  a km,  and  92  m for  a range 
of  10  km.  Neither  of  these  ranges  would  yield  much  of  a change  in  mean  sound 
speed,  there  being  only  a 0.1  m/sec  change  from  the  surface  value  in  the  latter 
case . 

Table  4.  Meteorological  variables  for  second  example.  (Actual  case  at  WSMR. ) 
Parameter  Surface  200  m-level 

Virtual  temp.  281.5  K 283.1  K 

East  wind  0 7 m/ sec 


South  wind 


0 


2 m/sec 
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However,  larger  gradients  are  possible.  Consider  the  case  described  by 
Table  4.  One  has  for  propagation  in  the  westerly  direction: 

c = 336.4  m/sec 
o 

c900  = 337.4  m/sec 
a = 1.49  x 10-’  m_1 
b = 10.4  x 10  ^ m ' 
d = 3.0  x 10  ^ m ' 

The  mean  speed  along  the  direction  of  propagation  is 

<v^>  = 336.4 [ 1 + 5.9  x 10_4  RJ ] , 

where  R is  the  range  in  km.  For  1 km,  ^v^>  = 336.6  m/sec,  and  for  10  km,  <v  > 

= 356.2  m/sec,  a 20  m/sec  difference  from  v . For  R = 1 km,  Z = 15  m,  and  for 
R = 10  km,  Z = 1.5  km  (higher  than  the  200  m level,  requiring  knowledge  of  the 
meteorological  parameters  above  200  m).  Since  most  temperature  inversions  do 
not  prevail  for  1.5  km,  the  R = 10  km  case  would  not  be  as  drastic  as  calculated 
above . 

This  method  of  approximation  to  the  correct  weighting  of  upper  layer 
meteorological  data  is  not  significantly  more  difficult  than  the  presently  used 
method.  For  a particular  ray,  if  a + b < 0,  then  surface  values  would  be  used. 
If  a + b > 0,  an  analysis  using  surface  values  could  be  used  to  determine  R,  and 
Eq.  (20)  would  yield  the  maximum  height  obtained.  The  meteorological  data  eval- 
uated at  one-third  this  height  would  then  be  used  for  hte  mean  sound  speed  de- 
termination . 

Refracted  propagation  also  opens  the  possibility  of  multiple  arrivals. 


Some  rays  leaving  the  source  at  smaller  elevation  angles  are  refracted  back  to 
the  surface  where  they  are  reflected  and  finally  arrive  at  the  receiver.  These 
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rays  arrive  later  than  the  main  refracted  ray,  and  with  amplitudes  dependent  on 
the  reflection  coefficient  at  the  surface  and  the  angle  of  incidence.  Since 
the  maximum  height  depends  on  the  square  of  the  distance,  the  correction  in  the 
time  of  arrival  depends  on  the  inverse  square  of  n + 1,  where  n is  the  number 
of  reflections.  Thus  the  travel  time,  which  for  the  main  refraction,  is  given 
by  Eq.  (B20). 


where  Z is  the  maximum  height  obtained  for  the  range  R. 

In  the  analysis  of  the  sound  ranging  problem  where  the  sound  speed  is  a 
function  of  range,  as  it  is  both  for  the  finite  amplitude  effect  and  the  re- 
fracted effect,  care  must  be  exercised  in  the  formation  of  the  equations  to 
assure  that  the  travel  times  reflect  the  mean  sound  speed  over  the  ray,  rather 
than  simple  time  differences  determined  by  the  local  sound  speed. 


h ns 
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III.  THE  RESPONSE  OF  THE  MICROPHONE 

The  T-23  microphone  is  basically  a pair  of  Helmholtz  resonators  with 
heated  wires  inserted  in  the  neck  connecting  the  two  resonators.  As  the  air  in 
this  neck  moves,  wires  change  their  temperature  due  to  convection.  These  wires 
are  part  of  a bridge  circuit  which  develops  an  unbalanced  electromotive  force 
due  to  the  change  in  resistance  of  the  wires.  Figure  7 indicates  the  physical 
layout  of  the  resonators. 

The  equivalent  electrical  circuit  for  the  double  resonator  is  shown  in 
Figure  8.  The  acoustical  inertance  of  a neck  is  labled  M,  the  acoustical  resis- 
tance of  a neck  is  labled  R,  and  the  acoustical  stiffness  of  a cavity  is  labled 
C.  It  is  apparent  from  Figure  8 that  the  bottom  resonator  acts  as  a series 
LC-tank  circuit,  with  the  top  resonator  acting  as  a low-pass  filter. 

From  measured  parameters  of  the  microphones  and  from  the  actual  response 
of  several  microphones  as  measured  in  a pistonphone,  the  frequency  response  of 
the  microphone  has  been  determined.  The  details  of  this  determination  are  con- 
tained in  Appendix  C.  The  resulting  response  is  illustrated  in  Figure  9.  It 
is  noted  that  the  response  is  sharp] v peaked  at  about  17  Hz,  with  a -6db/octave 
roll-off  for  low  frequencies  and  a -18db/octave  roll-off  for  high  frequencies. 
Comparing  this  response  with  the  spectrum  of  the  artillery  pulse,  as  developed 
in  Section  I,  confirms  that  the  microphone's  frequency  response  does  not  pass 
enough  high  frequency  information  to  faithfully  reproduce  the  signal. 

The  problem  faced  in  increasing  the  frequency  response  of  the  microphone 

*'****  ' ■’  •’  • • . . • J - • •«  r 

is  that  this  will  also  increase  the  noi  e in  the  recorded  signal.  Whether  the 
signal -to-noise  ratio  is  improved  by  an  extended  frequency  response  depends  on 
the  nature  of  the  noise  spectrum.  The  noise  spectrum,  of  course,  depends  on 
a lot  of  variables,  including  terrain,  wind  fluctuations,  and  man-made  noise. 


I 

figure  H.  Kquivalent  electrical  circuit  for  Hot-wire  microphone. 

M md  R are  the  acoustical  inertance  and  resistance  for  the  necks  and 
C is  the  acoustical  stiffness  of  the  cavities. 
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Figure  9.  Frequency  response  of  the  T-23  microphone. 
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Some  studies  of  low-frequency  noise  due  to  turbulence  indicate  that  the  spec- 
trum falls  off  on  the  order  of  6db/octave^’^.  On  the  other  hand,  man-made  noise 
due  to  battle  would  have  a different  spectrum.  Since  the  explosive  source  as 
considered  in  Section  I has  a period  approximately  in  proportion  to  the  diameter 
of  the  bore,  small  arms  fire  should  have  a spectrum  which  is  flat  to  about 
900  Hz  (.30  caliber). 

The  transient  response  of  the  microphone  is  more  illuminating  to  the 
problem  of  determining  time  of  arrival  than  the  frequency.  The  response  of 
the  T-23  microphone  to  a triangular  pulse  which  is  an  approximation  to  the 
explosive  source  shown  in  Figure  2 is  obtained  in  Appendix  C by  the  method  of 
laplace  transforms.  The  results  for  various  periods  are  shown  in  Figure  10. 

The  amplitudes  are  normalized  so  that  a pulse  length  of  1000  msec  (almost  a 
Heaviside  step  for  this  microphone)  will  yield  an  output  of  unity. 

As  shown  in  Section  I,  the  periods  of  interest  in  sound  ranging  vary 
from  10  msec  (light  artillery)  to  10  msec  (heavy  artillery).  As  seen  from  Figure 
10,  these  pulses  have  responses  varying  from  35%  to  70%  of  the  response  to  a 
1000  msec  pulse  of  the  -same  amplitude.  There  is  also  variation  in  the  amplitude 
of  the  negative  over-shoot  of  the  response.  These  variations  are  shown  in 
Figure  11,  where  the  amplitudes  of  the  positive-going,  and  negative  going 
responses  to  pulses  of  various  periods.  Also  shown  in  this  Figure  are  the 
approximate  periods  for  light  (105  mm),  medium  (155  mm),  and  heavy  (8  in)  artil- 
lery pulses,  lint  notes  that  the  valleys  have  somewhat  less  amplitudes  than  the 
peaks  and  that,  for  light  artillery,  besides  the  smaller  pressures,  the  response 
is  only  about  one-half  for  that  of  heavy  artillery.  This  is  due  to  the  restricted 
frequency  response  of  the  microphone. 


Period  (msec) 


Figure  11.  Response  to  unit  triangular  pulse  of  period  . The  solid 
curve  is  the  positive  peak  response,  and  the  dashed  curve  is  the  nega- 
tive peak  response  which  follows.  L,  M,  and  H refer  to  light,  medium 
and  heavy  artillery. 
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Of  course,  one  must  see  the  signal  in  order  to  determine  the  location  of 
the  source,  however,  the  time  response  to  the  pulse  is  critical  in  determining 
the  arrival  time.  There  are  at  present  four  points  on  the  response  curve  which 
are  used  to  determine  the  arrival  time: 

(1)  The  Break 

(2)  The  Pressure  Maximum 

(3)  The  Cross-over 

(4)  The  Pressure  Minimum 

The  orthodox  order  of  preference  of  these  points  is  (1),  (3),  (2),  and  (4).^ 

One  observes  from  Figure  10  that  the  break  is  difficult  to  determine  in 
that,  for  all  periods,  there  is  very  little  signal  for  the  first  three  msec. 

This  is  due  to  the  sharp  cut-off  on  the  high-frequency  side  of  the  frequency 
response  of  the  microphone.  Assuming  a noise  of  ±0.0 5 (S/N  = 26db)  added  to 
the  response,  the  approximate  deviations  in  determining  the  four  points  in  time 
may  be  found  to  be: 

Break  = ±2  msec 
Maximum  = 18  msec 
Cross-over  = ±3  msec 
Minimum  = ±10  msec 

confirming  the  orthodox  order,  although  the  cross-over  is  almost  equivalent  to 
the  break. 

Unlike  the  detectable  break,  which  lags  about  3 msec  behind  the  actual 
signal,  the  peak,  cross-over,  and  valley  lags  depend  on  the  period  of  the  signal. 
This  is  seen  in  Figure  12,  where  the  lag  of  these  phases  are  shown  as  a function 
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of  the  period  of  the  pulse.  It  is  seen,  however,  that  the  lags  are  fairly 
insensitive  to  the  period;  a change  in  period  of  20  msec  resulting  in  a change 
in  lag  of  the  valley  (total  length  of  signal)  of  only  10  msec.  This  insensi- 
tivity has  two  results:  first,  it  makes  it  difficult  to  determine  the  period 
of  the  pulse  and  thus  the  type  of  artillery  from  the  shape  of  the  signal. 

Second,  it  allows  the  use  of  any  of  the  three  alternate  phases  (maximum,  cross- 
over, and  minimum)  to  be  used  as  arrival  times  (assuming  the  same  phase  is  used 
across  the  whole  array).  The  maximum  differential  in  pulse  period  across  the 
array  for  a particular  pulse  would  be  for  heavy  artillery  at  short  range  (1  km). 
Table  2 shows  this  to  be  about  3 msec  across  the  entire  array  for  an  extreme 
flanking  angle.  This  results  in  only  a variation  of  1 1/2  msec  for  the  valleys, 
with  less  for  the  peaks  and  cross-overs. 

Thus,  the  use  of  a sharply  peaked  microphone  has  both  advantages  and 
drawbacks.  Nearly  all  signals  look  about  the  same,  which  allows  phase-compari- 
son of  signals,  but  the  differences  in  input  pulses  are  masked.  In  particular, 
the  initial  response  to  the  signal  is  very  slow,  but  the  lag  time  is  almost  in- 
dependent of  the  period  of  the  input  pulse.  All  this  reminds  one  of  Heisenberg's 
uncertainty  principle  — and  for  good  reason  — the  product  of  the  bandwidth  of 
the  microphone  and  the  uncertainty  of  the  arrival  time  of  a pulse  must  be  a 
constant  for  a particular  signal -to -noise  level. 
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IV.  RECOMMENDATIONS 

Each  of  the  foregoing  sections  are  theoretical  in  nature  and  suggest 
experiments  that  would  test  the  theory.  The  following  recommendations  stem 
directly  from  this  investigation: 

(1)  A STUDY  OF  THE  NATURE  OF  THE  SOURCE  SIGNAL.  Although  the  triangular 
shaped  waves  due  to  explosions  have  been  studied  for  sources  that 
range  from  very  large  nuclear  blasts  to  conventional  TNT  explosives, 
there  may  be  questions  still  unanswered  dealing  with  artillery. 

How  symmetrical  is  the  wave?  Are  there  secondary  shocks,  creating 
"N"  shaped  waves?  What  is  the  exact  relation  between  charge  and 
bore  to  the  period?  These  measurements  require  microphones  with 
wide  frequency  responses,  and  these  should  be  placed  from  about 
500  m to  several  km  from  the  sources.  Also  required  would  be 
meteorological  data  to  confirm  whether  the  propagation  is  diffracted 
or  refracted.  Some  of  the  information  sought  in  this  study  may  al- 
ready be  available,  but  unknown  to  the  author.  The  rise  times 
associated  with  the  received  signals  could  answer  the  question  of 
what  is  the  effective  specific  acoustic  impedance  of  certain  types 
of  terrain. 


(2)  A STUDY  OF  REFRACTIVE  CORRECTIONS  TO  SOUND  SPEED.  The  present 

method  of  weighting  the  vertical  meteorological  data  appears  aci  hoc . 
Section  II  contains  theory  which  indicates  a method  of  weighting  which 
could  be  programmed  and  the  resulting  locations  be  compared  to  the 
present  method.  Besides  this  test  on  already  existant  data,  study 


(1)  should  yield  actual  travel  times  if  timing  is  included.  This 
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timing  data  could  also  indicate  whether  the  surface  meteorological 
or  the  presently  used  200-m  data  is  the  more  proper  data  to  use  in 
analysis . 

(3)  NOISE  STUDIES.  Arrival  times  could  be  improved  by  extending  the 
response  of  the  microphone  and  by  increasing  the  dynamic  range  of 
the  recorder.  This  is  only  true  if  the  signal-to-noise  ratio  is 
not  significantly  decreased.  There  is  very  little  known  about  the 
noise  present  underneath  the  wind  screen.  Here  again,  a wide-fre- 
quency  response  microphone  is  needed,  along  with  wind  speed  data. 

As  mentioned  in  the  introduction,  an  investigation  of  semi-automatic 
determination  of  arrival  times  has  been  made.  This  was  done  with  an  HP  9825 
calculator  with  the  input  data  digitized  at  a 5 msec  sampling  rate  by  hand.  The 
calculator  was  able  to  detect  the  presence  of  signals  by  their  magnitude  to 
within  about  10  msec.  Closer  analysis  determined  that  although  the  breaks  were 
difficult  to  detect,  the  peaks,  cross-overs,  and  valleys  were  detected  to  within 
about  a msec.  In  addition,  cross-correlation  of  entire  signals  resulted  in 
arrival  time  differences  of  less  than  one  msec  except  for  difficult  signals. 

A system  compatible  with  such  a calculator  would  acquire  data  with  micro- 
phones with  a flatter  frequency  response,  and  thus  more  information  content. 
Since  the  algorithms  used  would  not  depend  on  hand  calculations,  the  microphones 
could  be  placed  in  non-linear  arrays,  enabling  distinction  between  forward  and 
rear  sources.  The  acquired  data  would  be  stored  digitally  on  magnetic  tape  and 
displayed  on  CRT's.  Operators  would  have  the  flexibility  of  filtering  and 
amplifying  the  data  for  display. 
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The  most  efficient  development  of  such  a system  would  be  by  a team  com- 
posed of  people  with  expertise  in  the  following  areas: 

(1)  Field  sound-ranging 

(2)  Microphones 

(3)  Computer  hardware 

(4)  Computer  software 

(5)  Atmospheric  propagation  of  sound 

(6)  Information  theory. 

Some  understanding  of  what  the  seismic  industry  has  accomplished,  and/or  famili- 
arity with  the  Navy's  SONAR  program  would  be  helpful  also. 

Finally,  it  is  quite  possible  that  the  biggest  problem  with  sound  ranging 
accuracy  is  the  fact  that  meteorological  data  is  obtained,  by  necessity,  out- 
side the  area  of  interest  and  previous  to  the  time  of  interest.  However,  this 
has  yet  to  be  proved,  and  therefore  improvement  in  time-of-arr ival  determination 
is  not  necessarily  an  academic  exercise. 
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APPENDIX  A 

THE  PROPAGATION  OF  A WEAK  TRIANGULAR  SHOCK 
Suppose  the  waveform  at  some  distance  r i from  the  source  is  similar  to 
Fig.  2 and  is  given  by 


P(x)  = < 


, x < 0 


P (1  - x/L  ) , 0 < x < L 

o o o 


, x < L . 

o 


(Al) 


As  the  shock  propagates,  the  pulse  length  increases  since  the  shock  at  x=0 
travels  faster  than  the  point  at  x = L.  The  rate  of  change  of  L is  given  by 


dL  = dL  dt  - I dL 
dr  dt  dr  c dt  ’ 


(A  2) 


where  c is  the  sound  speed.  The  velocity  of  the  shock  is  given  by 


O T ..+1  p 


(A3) 


where  y is  the  ratio  of  specific  heats,  P is  the  differential  shock  pressure 

and  p is  the  ambient  pressure.  For  weak  shocks  P/p  <<  1 and 
ro  o 


thus 


v - c(l  +^f> 

4y  P0 


dL  _ Y+l  P 

IF  - C(1  +T7r)  ■ c’ 

O 


(A4  ) 


(A5) 


since  the  end  of  the  pulse,  having  no  overpressure,  travels  at  the  speed  c.  Eq. 
(A5)  substituted  into  Eq.  (A2 ) yields 


dL  = _j+l  _P_ 
dr  4y  p 


(A6) 


The  total  energy  of  the  shock  wave  may  be  obtained  by  integrating  the 
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energy  density 


u - 2 , 

pc 

where  p is  the  density  of  the  medium,  over  the  entire  volume  of  the  wave. 
Assuming  spherical  propagation,  this  becomes 


(A7  ) 


E = 


/x 

a - ^2d*, 


2 2 

4;rr  p L 
3pc2 


(A8) 


(A9) 


If  energy  were  conserved,  one  could  solve  Eq . (A9)  for  P(r,L),  substitute 
into  Eq.  (A6)  and  integrate.  There  is,  however,  an  increase  in  entropy  across 
the  shock  which  depends  on  the  third  power  of  the  pressure  difference  across 
the  shock^ 


.0=2.  3_V  £ 

T 2 12  ’ 

\9P 


(A10) 


where  S is  entropy,  T is  absolute  temperature,  and  V is  volume.  For  an  ideal 
gas, 

/‘~v\  ..4-1  V 

(All) 


= Y+l  V_ 


o 2 
Y 

where  the  adiabatic  relation  pV ' = coust  has  been  used.  Substitution  of  Eq. 
(All)  into  (A10)  results  in 

3 


AS 


= nR(y+l)  /P 


12y 


03' 


(A12) 


where  the  equation  of  state  pV  = nRT  has  been  called  upon.  The  number  of  moles 


is  n and  R is  the  gas  constant. 
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The  rate  change  of  entropy  with  distance  is 


dS  = dS  ht  = l dS 
dr  dt  dr  c dt 


/ p \ ’ 

_ R(y+1)  I ] dn 

„ 2„  \p  / dt  ’ 

1 o,,  n • 


(A13) 


/P 

\ 

12y_C  " o 

where  dn/dt  is  the  rate  of  flow  (in  moles)  across  the  shock.  Again  assuming 
spherical  propagation, 

(A14 ) 


dn  2 cp 

dT  = A7Tr  “m  * 


where  cp  is  the  max  flux  and  M is  the  molecular  mass.  Thus 


dS  _ 4irr2(Y+l)P3 


dr 


2 9 

12y  P„  T 


where  the  ideal  gas  relation  p = pM/RT  has  been  used. 
The  rate  of  energy  loss  is  dE  = -TdS,  or 

,3 


dE 

dr 


47tr‘'(Y+l)  P~ 


12y2P  2 
o 


Division  by  Eq.  (A9)  leads  to 


1_  dE  = (y+1)pc~  P 

E dr 


2 2 

4YPo  L 


I oE 
E dr 


~ (— ) f , 

4 p L 
o 


where  c = y P /p  has  been  used  to  simplify.  Thus 


1_  dE  _ _ 1 dL 
E dr  L dr  ’ 


which,  upon  integration  yields  the  remarkably  simple  result 


(A15) 


(A16) 


(A17) 


(A18) 


(A19) 


J 
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EL  = constant. 


Thus 


2 2 2 

47Ir  “P  “L 

OOP 

9 

3pc 


2 2 2 
4nr  P L 


9 

3pc“ 


which,  when  substituted  int  Eq.  (A6)  becomes 


(A20) 


(A21) 


(A22 ) 


L 

o 


LdL 


. , P r 

Y+]  o o 

4y  p 

o 


(A23 ) 


L = { L [ L + 
o o 


Y+l 

2y 


r e (r/r  )]} 
on  o 


(A24  ) 


The  above  derivation  has  considered  the  increase  in  entropy  due  to  the 
finite  amplitude  of  the  shock,  but  has  neglected  the  increase  in  entropy  due  to 
viscothermal  and  molecular  relaxation  effects.  The  absorption  of  sound  in  air 
is  frequency  dependent  and  highly  sensitive  to  humidity.  High  frequencies  are 
attenuated  more  rapidly  than  low,  resulting  in  a rounding  off  of  the  initial 
pressure  peak. 
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APPENDIX  B 

MEAN  SOUND  SPEED  FOR  WAVES  REFRACTED  BY  CONSTANT 
GRADIENTS  IN  SOUND  SPEED  AND  WIND  VECTOR 
If  0 is  the  elevation  angle  of  a sound  ray  (normal  to  a wave  front)  and 
<f>  is  the  azimuth  angle,  then  the  law  of  refraction  yields^ 


= V 


where  0 and  $ are  the  original  source  values  for  the  ray  and  c is  the  phase 
velocity  of  sound.  The  first  equation  demands  that  the  azimuth  angle  is  constan 
and  the  second  is  Snell's  law  of  refraction.  In  a moving  medium,  the  phase 


velocity  is 


c = c + n-w  , 
P 


where  c is  the  sound  speed  in  a nonmoving  medium,  n is  the  unit  vector  in  the 


direction  of  the  ray,  and  w is  the  wind  velocity  vector  (the  velocity  of  the 
moving  medium).  This  is  in  contrast  to  the  group  velocity 

c = cn  + w , (B3) 

5 

which  is  the  velocity  of  energy  propagation. 

Since,  for  a particular  ray,  4>  is  fixed,  the  wind  vector  may  be  resolved 
into  two  components:  u in  the  direction  of  and  v normal  to  n.  Thus  Eq.  (Bl) 
becomes 


c 4-  u cos0 
cosO 


c + u COS0 
o o o 


If  the  atmosphere  is  considered  to  have  linear  gradients  for  c and  w, 
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c = c (1  + az) , 
o 


u = u + c bz 
o o 


v = v + c dz 
o o 


where  a,  b,  and  d are  constants  with  units  of  (length)  . Thus  Eq.  (B4)  has 
the  solution 


COS0  = COS0 


1 + az 

o 1 - bz  cos0 


If  the  additional  assumption  bz  cos0  <<  1 is  made,  then  Eq.  (B  ) becomes,  to 

o 

the  first  order 

cos9  = cos0  [1  + (a  + b cos©  )zl.  (B7) 

o o 

The  above  assumption  requires  that  the  total  change  in  wind  speed  over  the  path 
of  the  ray  is  much  less  than  the  sound  speed  itself.  Eq.  (B7)  may  be  solved  for  ; 


COS0  - COS0 
o 

(a  + b cos©  )cos0  ’ 
o o 


and  thus 


dz  = - 


sin6  do 

(a  + b cosO  ) cos0 
o o 


The  velocity  of  the  ray's  tip  is  given  by  the  group  velocity 

dx  _ a a 

— = c cos6  + u, 
dt 


(BIO) 


= c sinO  , 
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= c cost)  + U. 


(Bin) 


dF  = c sin0> 


where  x is  in  the  ^-direct ion , y is  normal  to  the  <J>-direction , and  z is  up.  It 
is  clear  that  the  ray  reaches  its  maximum  height  Z when  0=0,  and  that  the  path 
is  symmetric  with  respect  to  this  position.  Thus,  since 


dx  _ cos9  u 

dz  sin9  c sin0  ’ 


X = / + 11  dz 

“I  sin0  c sin0  ’ 


(Bll) 


(B12  ) 


where  X is  the  x-coordinate  of  the  point  of  arrival  of  the  ray  back  at  the  sur- 
face. (See  Fig.  Bl).  Eq.  (B9)  simplifies  this  to 


^ (a  + b cose  )cos0 
o o 


cosOdG  + I u 


(Bl  3) 


The  first  integral  in  Eq.  (B13)  is  simple  and  the  second  is  much  smaller 
(since  u/c<<1),  allowing  for  the  approximation 


u + c bz  u 
u o o ; o . . 

- = j t— r r — + bz. 

c c (1  + az)  c 
o o 


(B14 ) 


Eq.  (B14)  combined  with  (Bll)  yields 


> u 0 , 

X * t — r r — r r—  [sinG  + — ^ — - - + — — — (tan0  - 0 )].  (B15) 

(a  + b cos9  )cosO  o c a + b cos0  o o 

o o o o 
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Assuming  that  0 « 1,  this  simplifies  to 


x = !iL  2a  + 7A-  e 

a+b  |_  6 (a  + b ) 


2+^l  , 

o cl 


(B16) 


where  terms  beyond  0 have  been  neglected  in  the  trigonometric  expansions. 
Following  the  same  arguments. 


20  r , 2 v -i 

v = 2.  d 0 + _o 

a+b  6 (a  + b)  o c_ 

O 


(B17) 


a+b  |_6  (a  + b) 

where  Y is  the  net  displacement  of  the  ray  in  the  y-direction. 

Thus  the  x and  y components  of  the  ray  have  been  found.  The  mean  speed 
is  this  displacement  divided  by  the  travel  time.  The  ' ime  of  travel  is  obtained 
from  integrating 


dt  = 


dz 


c sin9  ’ 


(B18) 


or,  with  the  use  of  Eq.  (B9), 


T = 


2 f ° d9 

(a  + b cos©  )cos0  I c (1  + az) 
o o 1 


where  T is  the  total  travel  time.  This  may  be  approximated  by 


T = 


c (a  + b cosQ  )cos0 
o o o 


/ 


(1  - az)  d0. 


(B19) 


(B20) 


which  becomes,  with  the  use  of  Eq.  (B8) 
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T = 


c (a  + b cosO  )cos6 
o o 


[a (tan9  - 0 ) q 

, 2 

a + b cosO  J ’ 

- o J 


or,  for  small  0 , 

’ o 


T _ 20O  fl  . a + 6b  2 

1 C (a  + b)  6 (a  + b)  o 

— 


(B21) 


(B22 ) 


L. 
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The  mean  speed  in  the  x-direction  is  thus 


y 9 “ u 

<v  > = — = c 1+  — + — , 

x T o 6 c 

o 


(B23) 


and  the  mean  speed  in  the  y-direction  is 


,,  v d0 

i o.o 

<v  > = — = c — + — — 
y T o c 6 

o 


(B24 ) 


Tliese  equations  are  the  mean  speeds  for  a ray  leaving  the  surface  at  an  eleva- 


tion 0 and  arriving  at  (X,Y)  at  time  T.  The  results  depend  on  a linear  gradi- 


ent in  c and  w which  is  not  too  steep  and  they  are  correct  to  the  second  power 


One  must  now  replace  0 ^ with  the  proper  function  of  the  distance  between 


source  and  receiver,  since  is  not  generally  known.  The  total  distance  tra- 


velled by  the  ray  (measured  along  the  surface)  is 


R = (X“  + Y 


„ . tv  r-  0 , ,,  0 u -i  ‘-i 

" a + b 6(a  + b)  o coJ  ’ 


(B25) 


correct  to  order  0 and  to  first  power  in  u /c  . Thus 
o o o 


0 = 4(a  + b)R, 

o 


(B26) 


correct  to  order  R*-.  This  yields 


<v  > = c [1  + — + -rr  (a  + b)  R ] , 
x o c 24 
o 


(B27 ) 


- r° 

- co 

n 


+ (a  + b)R2 


(B28 ) 


The  maximum  height  obtained  by  the  ray  is  given  by 


1 - COS0 


(a  + b cosO  )cos0  * 
o o 


(B29) 
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} 


or,  for  small  0 , 


6 


2 (a  + b) 


(B30) 


Thus  retaining  terms  of  order  0 ~/(a+b),  which  has  been  done,  corresponds  to 
retaining  terms  linear  in  z,  consistant  with  the  constant-gradient  atmospheric 
model.  Finally,  it  may  be  noted  that  Eq.  (B30)  may  be  written 


Z = I (a  + t)R2,  (B31 ) 

giving  the  maximum  height  in  terms  of  the  range  and  the  gradients. 

The  substitution  of  Eq.  (B31 ) into  (B27)  yields 

u . 

<v  > = c [1  + — + ^ (a  + b)Z], 
x o c 3 
o 

= c(|)  + u(f),  (B32) 

or,  the  mean  speed  in  the  x-direction  is  the  speed  evaluated  at  Z/3,  one-third 
of  the  maximum  height.  Likewise  Eq.  (B28)  becomes 

v 

<v  > = c (—  + — dZ] 
y o c 3 
o 

= v (~)  (B33) 

or,  the  cross-wind  evaluated  at  1/3  of  the  maximum  height. 
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APPENDIX  C.  T-23  MICROPHONE  FREQUENCY  AND  PULSE  RESPONSE 

The  T-23  double  resonator  hot  wire  microphone  has  the  acoustical  configu- 
ration shown  in  Figure  7.  The  analogous  electrical  circuit  is  shown  in  Figure  8. 
This  circuit  is  composed  of  elements  defined  to  be: 


M = Inertance  of  neck, 

R = Acoustical  resistance  of  neck 
C = Compliance  of  cavity, 

= V/pc“  = V/yp 


where  y is  the  ratio  of  specific  heats  and  p is  the  ambient  pressure.  The  sub- 
scripts refer  to  the  upper  resonator  (1)  and  the  lower  resonator  (2).  The 
acoustic  pressure  incident  on  the  microphone  is  P(t)  and  the  volume  velocity 
flowing  through  the  necks  is  U(t). 

For  P(t)  = P^e1'*",  the  loop  equations  are: 


P = iuiM.U.  + R.  U 
o 11  11 

0 = iuM.,U0  + R0U0 


U1  -U2 
u)C, 


1 


(Cl) 


where  and  U0  are  the  amplitudes  of  the  volume  velocities  through  the  necks. 

Since  the  output  is  proportional  to  U0,  the  relative  response  of  the 
microphone  is  given  by 


1 

Z ’ 


(C2 ) 


where  Z is  the  effective  acoustical  impedance.  The  solution  to  the  set  of  equa- 


tions  (Cl)  yields 


where 


Z = R'  + iu)M'  - i 


R'  = Rl(1  + - u2M2C1)  + R2(l  - u,2!^). 


c. 

M'  = Mjd  + " 2M2Cl)  + M2  + R1R?C1. 


The  output  voltage  is  thus. 


E = 20  log (K/ | Z | ) db. 


where  K is  a constant  of  proportionality  and  jz|  is  the  absolute  value  of  the 
complex  impedance.  The  phase  angle  of  the  output  with  respect  to  the  input  is. 


where , 


0 = tan  (-X'/R'), 


X'  = uiM’ 


The  acoustical  impedance  of  a small  tube  is 


, L / 8u  , 4 \ 

Zt  = — (~  + 3 la)P  • 
irr  'r  / 


where  L and  r are  the  length  and  radius  of  the  tube,  P is  the  viscosity  of  the 
air,  and  p is  the  density  of  the  air.  The  criteria  for  smallness  depends  on  the 
thickness  of  the  viscous  boundary  layer  which,  in  turn,  depends  on  the  frequency 
of  the  signal.  In  addition,  the  effect  of  the  impedance  miss-match  between  the 
infinite  medium  and  the  tube  is  sometimes  approximated  by  adding  a correction  to 


the  length  of 


(CIO) 
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The  volumes  of  the  cavities  and  the  characteristics  of  the  necks  have  been 
measured.-  Using  C = V/Yp  and  Eqs.  (C9)  and  (CIO), 

C^  = 1.9  x 10  3 cm^/dvne 

= 0.7  x 10  3 cm3/dvne 

R^  = 3.4  acoustical  ohms, 

R?  = 30  acoustical  ohms, 

= 3.8  x 10  - gm/cm^ 

M,  = 8.2  x 10  gm/ctn* 

In  the  above,  the  characteristics  of  the  25  Hz  plug  in  the  top  neck  has  been 
assumed.  Since  the  tubes  in  this  case  are  somewhat  thicker  than  the  boundary 
layer,  the  values  for  R and  M should  be  reduced. 

3 

Figure  Cl  shows  the  relative  response  of  several  T-23  microphones. 

This  response  was  measured  in  a pistonphone  and  is  accurate  to  within  about 
±1  db.  The  solid  line  beyond  about  17  Hz  is  calculated  from  the  above  values, 
and  is  seen  to  be  an  accurate  representation  of  the  measured  response.  For  less 
than  about  17  Hz,  there  are  two  lines  shown.  The  dashed  line  is  the  calculated 
response  using  the  values  determined  above.  The  solid  line  is  an  empherical  fit 
to  the  data  using: 

Rj  = 2.4  acoustical  ohms, 

R0  = 20  acoustical  ohms, 

= 3.6  x 10  £ gm/cm^ 
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M,  = 7.8  x 10  " gm/cm  , 

assuming  that  R and  M should  be  decreased  due  to  the  finite  radius  of  the  tubes. 

As  may  be  seen  from  Figure  Cl,  the  above  empirical  values  closely  fit  the 
observed  response  curve.  However,  in  the  case  of  sounding  ranging,  the  impor- 
tant response  is  charactei ized  by  the  response  to  a triangular-shaped  explosive 
pulse,  rather  than  a sine  wave.  The  response  to  a pulse  is  obtained  from  a 
Laplace  transform  of  the  circuit.  This  is  obtained  from  the  frequency  response 
by  substituting  s for  iw: 

rjffen  • »3<",m2ci>  + »2fV2ci  * Wd  + sWi  * "iVc2 + rir2ci+m2) 


+ R1  + ElVC2  + sfc'' 


(Cll) 


where  L[  ] indicates  the  Laplace  transform.  Some  algebraic  manipulations  yield 


2 2 

C w.  cu  s 

[E(t>]  = L[P(t)]. 


(C12 ) 


where  ^ ^ . and 


f(  ) = s*  + (fj  + 5,,)s^  + (ti^“  + u>0“  +u>|)"  + 


+ (S.uij  “ + j ) s + Uj 


(C13) 


w 2 = l/M2Cr  ^ = R]/M1  and  62  = R^L,. 

If  one  defines  x = s/ui  ^ » then  Eq.  (C13)  becomes 

U U 3 2 

g(x)  = o^fx4  + a^x  + a^x  + a^  + aQ], 


(C14) 
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where 


We  wish  to  factor  the  quartic  Eq.  (C14)  by  Ferrari 's  method, 
to  obtain  the  resolvent  cubic 


3 2 

y - a2y  + (a^ 


2 2 

4a  )y  - (a,  “ + a a 
o 1 o 3 


4a  a 
o 


2 


For  the  case  at  hand, 

a3  = 4.09, 
a9  * 8. 54, 

aj  » 10.12, 

a = 5.88, 


which,  by  numerical  methods  yields  a real  solution  y = 4.857. 
results  in  resolving  the  quartic  into  two  quadratic  equations 


w 


12 


(xJ  + Cj  (+)x  + c2(+))(x2  + Cl(_)x  + 


(C15) 


The  first  step  is 


(016) 


This  solution 


2(_)),  (Cl  7) 


More  algebraic  manipulations  result  in 


e (s) 


sp(s) 

^ (s)4>0  (s) 


for  Eq.  (C12 ) where  e(s)  = L[E(t)],  p(s)  = L[P(t)], 

♦ j = (x  + a1 )“  + b1 
4>2  = (s  + a.,)“  + b.,~ , 

and  where 

a1  = 112.9, 
a 2 = 55.0, 

bj  = 52.5, 

b,  = 120.1. 

Assuming  a triangular  pulse  for  P(t), 

fo  , t < 0 


and 


P(r) 


1 - t/T  , 0 > t > T 


I 
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then 


£[P(t)] 


o 


t/T)e  Stdt, 


or, upon  integration. 


P(s)  = 7 (1 


'S 


-s  T 


). 


(C20) 


(C21) 


Substitution  into  Eq . (C18)  yields 


r- 


e(s)  = (s)  + e.,(s)  + e^(s),  (C22) 

where 

el(s)  = i>]  (sH,,(s)  ’ 


e2(S>  -'s4>1  (s)i>2  (s)  ’ 


e~sT 

e3(s)  = "^oTH.7(7y  • 


Considering  e^Cs)  first,  this  may  be  written 

A + B s A + B s 

ei(s)  =-|^7)-  + -7^)- 


(C23) 


A table  of  Laplace  transforms  then  yields 


-v 


E1(t)  = 


blCl 


2 tPrlsin  bLt  - Pncos  b^] 


-a2t 


+ ^ [Pr2sin  b2c  “ P12COS  bTt^’ 


b2C2 


(C24 ) 


whore 


9 9 9 


>Vl  = (a2  ■ V + V - bl  ’ 


Pil  = 2bl(a2  _ V’ 


9 9 9 


(a2  “ ai)  + b2*" + bi  ■ 4brb2* 


9 9 9 


Pr2  = (al  ' a2r  + b2~  + br  ' 4bl*h2"  ’ 


Pi2  = 2b2(ai  " a2>» 


c2  = C1 


This  portion  of  the  response  of  the  microphone  is  that  due  to  a unit  step  function, 
and  is  shown,  normalized  to  unity,  in  Fig.  C2. 

Next,  considering  e~(s). 


e (s)  = - .7-  e.  (s) . 

2 .si 


(C2-1) 


which  becomes 


E.,(t)  = - - 


K j (y)dy. 


(C26) 


upon  taking  the  transform.  Integration  yields 


Ft  (t ) 


- 2 I(Bbl  _ Aa1)sinb1t 

lar+br 

Ab  + Ba1 

- (Ba j + Ab  )cos  h^t)  + — ^ j 

ar + v 


(Eq.  cont.  on  next  page) 


(C27) 


where,  for  this  case. 


a9  + b 9 


— [ (Bb.,  - Ca  , )sin  b^t 


Cb,  + Ba 


- (Ba?  + Cb^cos  b,t]  + — j =•  V , 


a2  + b2 


A = 6.19 
B = 2.51 
C = 1.50. 


This  function,  the  response  to  the  triangular  portion  of  the  pulse,  is  shown  in 
Fig.  C3,  with  the  same  normalization  as  Kj(t). 

Finally,  the  last  term  may  be  resolved  bv  the  use  of 


F(t 


') 


. 1 — . s 
L e 


f(s). 


(C28) 


yielding 


E3(t)  = 


t ^ T 


-F-2  (t  - T)  , t > T 


The  required  function  E^t)  and  F.o(t)  are  tabulated  in  Table  Cl. 


(C29) 
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Table  Cl. 

Response  functions 

so  that  E,  (T) 

1 max 

duration  T is  E(t) 

t < 7,  and  E (t)  = 

for  the  T-23  microphone. 

1.000.  The  response  to 
= E (t)  + E2(t)  + E3(t), 
-E_  (t  - '.")  for  t > T. 

Functions  are  normalized 

a triangular  pulse  of 
where  E^ (t ) = 0 for 

EL(t) 

7E0(t) 

tlra  (msec) 

(volts) 

(volt • sec ) 

0 

.000 

.0000 

5 

.076 

-.0001 

10 

.376 

-.0012 

15 

.74 

-.0040 

20 

.97 

-.0084 

25 

.97 

-.0133 

30 

.77 

-.0177 

35 

.47 

-.0208 

40 

.18 

-.0224 

45 

-.02 

-.0228 

50 

-.  12 

-.0224 

55 

-.14 

-.0217 

60 

-.  10 

-.0211 

65 

-.05 

-.0207 

70 

.00 

-.0206 

75 

.03 

-.0207 

80 

.03 

-.0208 

85 

.03 

-.0210 

90 

.01 

-.0211 
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